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Abstract 



A solution of the linearized Einstein's equations for a spherically symmetric perturba- 
tion of the ultrarelativistic fluid in the homogeneous and isotropic universe is obtained. 
, Conditions on the boundary of the perturbation are discussed. The examples of particle- 

00 ' like and wave-like solutions are given. 

O 

O 1 Introduction 

\o ' 

The description of metric perturbations in relativistic cosmology is significance in the theory 
of creation of the universe large-scale structure. Spherically symmetric perturbations are an 
important and interesting class of that perturbations. 
bX)' It is easy to understand that the spherically symmetric perturbation may be represent as 

outgoing and ingoing waves (traveling from and into the center of the configuration) . The ingo- 
^ ! ing wave amplitude increases as soon as the wave approaches to the center of the configuration. 
^ I Therefore the nonlinear stage of the development of the spherically symmetric perturbation 
must appear more quickly in the contrast to the case considered by Lifschitz [|l| for flatwave 
perturbations. As the wave velosity in the ultrarelativistic fluid is equal to the velosity of sound 
c/ so the time T of appearing of the nonlinear stage is only determined by the typical size 
L of an initial perturbation: T = \/3Lc~^. Below we describe this process explicitly. 

2 Ultrarelativistic fluid in spherically symmetric space- 
time. Rigorous equations 

The spherically symmetric metric may be taken in the form 

ds^ = e^dT]^ - [dr^ + (dO^ + sin ^6 d^^)] , (1) 
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where the metric parameters v and A are the functions of the time coordinate r] and a radial 
space coordinate r . With this notations the nontrivial Einstein's equations take the form 

(0 = ^ = l,£ = 3p) 



^A'2 + AV + -(A' + i/') 
2 r 



2(A" + i.") + i^" + ^(A' + J^') 
1 

-A / \// 



.. 1 . . 3.2 

A--A.+-A 



87rp (l + 4i;2) , 



3 .2 



— A -- Az> +7 A = 87rp, 



(a" + 1a- + ^A')+^A^=8.p(3 + V), 
(2 A -v' a) = 87re^''-^)/^4pi;VrTi^, 



2x/3 



(2) 

(3) 
(4) 
(5) 



where e is the energy density and p is the pressure; v — vFe^l'^ is the frame projection of radial 
velisity and vJ' is the radial coordinate of 4-velosity of fluid; a dot and a prime denote partial 
derivatives with respect to 

and r, respectively. 

There are two consequences of these equations 
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_A" - v" + ^A' ' + AV - ii/' 2 + 1 (A' + v') 

Zi Z f 



Sttv 4p, 



(7) 



— e 



2A" + u" + -(X'^ + AV + u'^) + - (2A' + u') 

+6-^^ {'\\\ -\ \v \ = 0. 



(8) 



3 Background spacetime 



In isotropic spherical coordinates the metric of the background space-time is 

1 



ds —a Idri j 



where 



y = Vl + Kr%-^, 



(9) 



(10) 
(11) 



aoT, for K — 0; 

a = < aosinh("\/l2T6~^), for K = -1; 

^ aosin(VT2T6-^), for X = +1, 

K = 0, ±1, b, oo are the constant parameters of the background space-time. 

The energy density £0, the pressure po and the frame projection of radial velosity of the 
backgound fluid Vq are 

,3 / o o \ 



£0 = 3po = ^ 
Stt 



4^( y"--y" + -y' ]+-,a 
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(12) 
(13) 
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4 Perturbation equations and general solution 

Let us define small perturbation quantities as follows 

5e = ?>5p = e — Sq, (14) 
5v = v-\n [a^] , (15) 

(16) 

We will assume that the frame projection of radial velosity of fiuid v is the small quantity 
too. Substituting relations (|l^p!6D into Eq.(|^) and keeping only the first-order quantities in 
perturbation we obtain 



5\ = \-\n (a^/y^) . 



{6\" + 6u") + - {6\' + 6u') - ^ {6\' + 6u') = 0. 
r y 



(17) 



The solution of this equation is 
5\ + 5v = 



Ci + C2r\ 

C3 + CJ^-'K-^y-^ 



for K = 
for K = ±1, 



where Ci, C2, C3 and C4 are the arbitrary functions of the time coordinate r. 
If we require that 

6\ = = for r > tq^t), 



(19) 



where vq is the some function of r, i.e. the space-time is asymptotically homogeneous and 
isotropic, then 

Ci = C2 = C3 = = (20) 

and 

5\ = -5v. (21) 
Substituting the Eqs.(|l5D, (|l|) and (|lD into Eq.(|), we get in the linear approximation 



-y 



5\" + 2(- - y^)5\' 
r y 



+ 5X + A-6X + A-6X = 0. 
a a 



(22) 



If we introduce new radial coordinate 

r 

dr 



L{r) 



l + K{r/hf 
and make the substitution 



6Artanh(r6 ^) 
6arctan(r6~^), 



K = 0; 
K = -l- 

K = +l 



ar OT \a/ 



where (p = (f{T, L), then the equation ( p2D can be rewritten in the form 



d_ 



(9V 4:K 
"dl? ~ "62 



(23) 



(24) 



(25) 
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One can easily integrate this equation 



9V 4ir 



-UrV = aF (L) 



(26) 



where F (L) is an arbitrary function of L. 

A general solution of Eq. (^) is the sum of the general solution of the equation 



9V 4ir 
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(27) 



and the particular solution of the equation (|26|) . But this particular solution does not change 
the form 6X, because it can be rewritten as aW, where W = W{L) and 



d'^W 16K 
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W = F. 



{21 



Therefore we can substitute the solution of equation (|27D but not equation (p6D into the ex- 
pression (pW). In the case K = this solution is 



^{t,L) = $+(r + L) + <l>_(r-L) 



(29) 



and in the case K = ±1 



T+L 



dt^ 



Joi-Jir-L) (r + L-t)]dt 



+ 



dt^ 
+ 6 



Jo 



{t + L){t-L- t)^ dt 



Jn 2t- 



J |t=o 
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(30) 



where $+, / and g are arbitrary functions and Jo is the Bessel function. 

Thus the relations (p^ and (^TD describe the spherically symmetric perturbation of the 
metric satisfying the condition (|T9D on the background of spacetime (|^), if we take into account 
the relations (H) and (|30|). 

The perturbations of the density 5e and the frame projection of the radial velosity v are 



a ■ a \ V 
87rfe = 87r3(5p = ^ U A - ^ 



\r y I 



Sttv 



4en 



+- [y -2 — + -y SX 

y \ y r I 



6X+-SX'). 



(31) 
(32) 



In the case 



K = 0, $+ (x) = D+ + Di x + D2 + Ds x^ $_ (x) = + Di x - D2 + x^ (33) 
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where D+, Di, D21 -D3 are constants, we obtain a Newtonian potential 

5\ = 

caused by a particle of variable mass 



^ 2m 

du = 

ar 



m = - ^{D+ + D_) + —2D3 

a 



2a? 



(34) 



(35) 



In the conclusion of this section we would like to note that the perturbations in the homo- 
geneous and isotropic universe are gauge-dependent 0]. Therefore, we should know whether 
the perturbative quantities 5e, v, S\, 81/ obtained in this section are the actual physical pertur- 
bations or merely an artifact of gauge. 

In concequence of the symmetry of problem the perturbations described in this paper are 
scalar ones. The gauge-invariant amplitude of density perturbation with the notations of the 
paper [ffl is 



5_3(l + w)l|fi?W--4°) 



kS 



(36) 



The choice of the isotropic coordinates (|^) in our paper corresponds to a longitudinal gauge. 
With the notations of the paper H this give 



5(0) = H. 



(0) 
T 



0. 



By comparing the corresponding expressions in two papers we find 

The relations between the other perturbative quantities can be determined analogously 

g(0) ^ ^(O)g(O) ^ 



k ^ 



(0) 



k k5 



1 ( frW . S • (0)' 



g(°) = AQ 



(0) 



5u 



g 



(0) 



(37) 
(38) 

(39) 

(40) 
(41) 



iFiom this relations we find that the perturbation quantities introduced in this section are 
actualy physical as well as the gauge-invariant quantities egQ^^\vPQ^^\ ^aQ^^\ ^hQ^^^- 



5 Boundary conditions 



The solution (34) is physically unacceptable for the description of the gravitational field of 
perturbation which arose as a result of fluctuation at the moment of "time" r = Tq, since it 
is inconsistent with the principle of causality. Therefore, the boundary conditions should be 
formulated in the way that at least beyond the light horizon, the potential SX should vanish 
together with its derivatives. Such boundary conditions are in according with the "birth" of 
a perturbation as a result of the redistribution of Robert son- Walker matter. In fact, however. 
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the horizon of the perturbation is not the hght cone but the sound cone, since the density 
perturbations extend with velosity of sound. Thus the boundary conditions at the sound 
horizon L = t — Tq are taken in the form 



L=r-To 



0. 



L=T-TQ 



Substituting Eqs. (p9D and (0) into this relations, we obtain in the case K = Q 

$'+ (2r - To) + (0) - - [$+ (2r - tq) + (0)] = 0, 

r 

$'1 (2r - To) - $ (0) - (2r - tq) + ^ [<f + (2r - tq) + $ _ (0)] = 0, 



r 



$'| (2r - To) + (0) - - [$V (2r - tq) + (0)1 + - [$+ (2r - tq) + (0)] = 



and in the case K = ±1 



(42) 



(43) 

(44) 
(45) 



(2r - ro) + (0) - - {/" (2r - To) - /" (0) + h [f" (0) + g'" (0)]} 



(46) 



(2r - ro) 
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[/'" (0) + g'" (0)] 



- g"" (0) - - [/'" (2r - ro) - g'^ 
a 

'2r2 a2r\ ^„ , , a ^, , , /2r a \ , , 



T-O)-/ 



0, 



(47) 



(2r - ro) + g"" (0) + + — j [/'" (0) + (0)] + ^g'" (0) 

+ 1 (1 + 6^) [/" (2r - ^o) - /" (0)] - (0) - (0) 

+ ^[/(2^-^o)-/(0)]=0, 

where a prime denotes a derivative of the function with respect to its argument. 



(4J 



6 Case of spatially flat universe {K = 0) 



Integrating Eqs.(|43|)-(P5|), we obtain 

$+ (x) = 



(ro) - X (ro) - ro$': (ro)] - (ro) + ro<l>'_ (ro) - $_ (ro) . (49) 
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We denote a new function 



^ (r - r - To) = $_ (r - r) - $_ (tq) - ^<^'_ (tq) ■ (r - r - tq) 



2! 



<f - (ro) ■ (r - r - tq)' 



Then it is possible to reduce the relation (^) to 



6X 



ar dr 



^ (r - r - To) 



(50) 



(51) 



(52) 



and the conditions (^31) -(^Sj) to 

^ (0) = ^' (0) = ^" (0) = 0. 

If the function \E' (x) satisfies to the condition 

^ (x) = for a; < 0, (53) 

then the relations (|5TD-(p2D describe the spherically symmetric fluctuation of gravitational field 
in region t — r — tq > that satisfies to the conditions 



6X\ 



d_ 

dr 



6X 



d_ 



6X 



0. 



(54) 



r = T — TO 

For more detailed study of the solution near the boundary we represent the function \E' (x) 
as a power series 

^(x) = 5]^„a;"if(x), (55) 



n=0 



where are constants and 



H(x) 



0, X < 0, 

1, X > 0. 



Taking the relation (|5^ into account, we obtain 

vl/g = v]/^ = vl>2 = 0. 

Then relation (^) reduces to 



1 ^ 



(5A = — V — (r-r - tq)" ^ ( n - 1 
The solution obtained in 

5X = ^3 



r + To 



(56) 



(57) 



(5^ 



2rV 1 3 r2 



+ 



a \ar 



2ar 2ar3 



H(t - r) 



(59) 



corresponds to the first term of the series (|58|) for tq = 0. The solution (^) is a particle-like 
one if the function 



m 



d_ 



^ (r - Tq) 

2a 



(60) 
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is not equal to zero. This function describes the mass of particle in the center of configuration 
and in the case considered in B it is 



m 



The wave traveling from the center of configuration 

^0 



SX 



arr 
0, 



47r 27r 1 2it 1 

— sm —X H — cos —X 

/ I r I T 



2n 

1 — cos —X 



for < X < I; 
for X < 0, X > /, 



(61) 



(62) 



where x = r — r — tq, \l/o and / are constants, may be an example of nonparticle-like solution in 
the case t > tq + I. The function \1/ that corresponds to that solution is 




for < X < /; 
for X < 0, X > /. 



(63) 



7 Conclusion 

We have considered the spherically symmetric perturbation of the cosmological perfect fluid 
with the equation of state p = e/3. The general solution of the linearized Einstein's equations 
which describe this perturbation has been obtained. The solution is represented as outgoing and 
ingoing waves moving with velosity of sound. In order to consider the perturbation in causally 
connected region we have imposed the boundary conditions on the solution. In this case we 
have investigated the behavior of solution near the boundary which is the sound horizon. We 
have given the examples of particle-like and wave-like solutions. 

We would like to finish our article by the following remark. The point of view that the 
large-scale universe structure is pancake-like is well-known |^. This point of view is based on 
the fact that fiat perturbations (with the wave length more than the Jeans one) begin to be 
increasing only with the decoupling moment on the nonrelativistic stage of the universe exten- 
sion. However, the obtained solution shows that spherical perturbations increase already on 
the ultrarelativistic universe extension stage. Thus it may be turned out that at the decoupling 
moment the spherical perturbations will have a bigger amplitude than the fiat ones, and in this 
case we will have to revise the scenario of creation of large-scale universe structure. 
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